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^ Abstract 
H 

We study thermal neutralino dark matter in an effective field theory extension of the 
MSSM, called "Beyond the MSSM" (BMSSM) in Dine, Seiberg and Thomas (2007). In 
this class of effective field theories, the field content of the MSSM is unchanged, but the 
little hierarchy problem is alleviated by allowing small corrections to the Higgs/higgsino 
part of the Lagrangian. We perform parameter scans and compute the dark matter relic 
density. The light Higgsino LSP scenario is modified the most; we find new regions of 
parameter space compared to the standard MSSM. This involves interesting interplay 
between the WMAP dark matter bounds and the LEP chargino bound. Wc also find 
some changes for gaugino LSPs, partly due to annihilation through a Higgs resonance, 
and partly due to coannihilation with light stops in models that are ruled in by the new 
effective terms. 
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1 Introduction 

In the next few years, there will be a significant amount of new data in dark matter 
physics and TeV scale particle physics. If low-energy supersymmetry is relevant to the 
real world, and in particular if it provides dark matter with the right properties, it is 
important to keep asking the question as to what extent the Minimal Supersymmetric 
Standard Model (MSSM) is the best framework for analyzing the new data. 

The MSSM is "minimal", in that it is the minimal supersymmetric completion of 
the standard model, with the only supersymmetry breaking coming from operators of 
dimension three or lower. From the low-energy point of view, this is a very special imple- 
mentation of low-energy supersymmetry. For example, the MSSM quartic Higgs terms are 
completely determined in terms of known electroweak gauge couplings by supersymme- 
try, and so are small (of the order A ~ 0.07). This is quite unlike the nonsupersymmetric 
standard model, where the quartic Higgs coupling A is fairly unrestricted, and need not 
be small. This is one reason why the nonsupersymmetric standard model has no problem 
with the present LEP lower bound on the lightest Higgs mass, whereas it is generally 
considered to be somewhat of a problem for the MSSM (more on this in the next section). 

Part of the reason for imposing this simplicity at low energy in the MSSM is standard 
coupling unification, where the "energy desert" (the assumption that essentially no new 
physics enters between the TeV scale and the GUT scale) means that all new operators 
due to new physics are suppressed by the enormous scale Mqut- Another reason to impose 
simplicity is purely practical — there are already over 100 free parameters in the MSSM, 
so why add even more uncertainty? 

Indeed, if the MSSM holds up to experiment, there is no need to add more uncertainty. 
However, what if it does not? For example, the recent fiurry of model building around the 
PAMELA results [IJ is partially motivated by the fact that the standard neutralino does 
not seem to be able to reproduce the PAMELA results. We do not attempt to model any 



3 



aspects of any particular new experiment in this paper, but the model building motivated 
by the new PAMELA data provides an example of how we may need to go beyond the 
"maximally minimal" framework for neutralino dark matter when attempting to ascribe 
any new phenomena to dark matter physics. 

But with guidance from experiment still scarce, how can we possibly generalize the 
MSSM, with its 124 parameters, in a meaningful way? The typical approach in particle 
phenomenology is to consider non-minimal extensions of the MSSM, where new fields and 
new parameters are added to the theory to solve specific problems, with some success 
and some degree of arbitrariness. Since so little is known about this non-minimal physics, 
it seems meaningful to attempt to be systematic about this generalization. A first step 
could be not to add any new particles at all ("minimally non- minimal" ) , but only allow 
indirect effects by any new heavy particles on the parameters in the Lagrangian, i.e. think 
of the MSSM as an effective field theory with a "scale of new physics" M <C Mqut, which 
can be as low as phenomenologically allowed]^ Even within this limited framework, there 
is already a bewildering variety of interactions and couplings one could add to or modify 
in the MSSM Lagrangian. 

A very rough way to classify such new parameters is by whether turning them on 
tends to make experimental constraints harder or easier to satisfy. For example, most 
parameters that already exist in the MSSM-124, such as off-diagonal squark masses, 
immediately bring the model beyond experimental bounds on flavor and CP violation 
when turned on to sizeable values. Therefore, one useful way to make the MSSM slightly 
non-minimal would be if we could identify some class of new parameters that rule in a 
given model, instead of ruling it out. Such parameters were identified by Dine, Seiberg, 
and Thomas [7]j^ They enumerated some operators at dimension four, five and six that 
go beyond the MSSM (hence the acronym BMSSM). These operators can be added with 
small coefficients to capture possible non-minimal physics at the scale M, which we take 
to be M > 5 TeV or somewhat higher. In particular, we impose that the scale M is 
sufficiently high — or the underlying microscopic model sufficiently restricted — that we 
can neglect operators suppressed as or more. By the logic above, we further restrict 

the operators at order 1/M by not turning on operators in the squark and slepton sectors 
that would worsen problems with flavor- changing neutral currents. This leaves only two 
operators in the Higgs sector, so two new parameters ei and €2. (It would be interesting 
to go beyond our analysis by performing a detailed study of flavor physics, CP violation 
or operators, but we leave that to future work. See also section^) 



^In the past one could have tried to argue against such new thresholds based on intuition from string 
models, but with the advent of stabilized string models with many interesting thresholds far below Mqut 
(e.g. [2J[3]), this seems overly restrictive. Incidentally, such string models often sport non- neutralino dark 
matter, as in |31|S], or general non-thermal dark matter, as in [5]. 

^For earlier work in specific models that identified similar operators, see [5J |5] and [TD] . 
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Roughly speaking, if the MSSM was aheady "natural", then adding some effective 
operators as small perturbations would not affect it much. It is a sign that the MSSM 
is very special that this small addition can change the way data would be interpreted 
in important ways. For example, tree-level deviations from the supersymmetric quartic 
Higgs coupling A could rule in models that were previously ruled out by LEP bounds. 
These models may have new phenomenological properties, e.g. for dark matter. Other 
aspects of the new models will be discussed in detail in this paper. 

The kind of dimension-four operator that shifts the quartic Higgs coupling has been 
analyzed before (e.g. [H]), but not together with the "companion" dimension-five non- 
renormalizable Higgs-Higgs-higgsino-higgsino interactions we need to include here. These 
operators are "companions" because they are both effective dimension five, i.e. the new 
quartic Higgs terms have coefficients that scale as 1/M. (There is of course another mass 
scale in the numerator to make the quartic Higgs coupling dimensionless, but this is a 
smaller scale intrinsic to the low-energy theory, such as the supersymmetric Higgs mass 
parameter /z.) The new quartic terms include a term that breaks supersymmetry and is 
technically hard. The new nonrenormalizable Higgs-Higgs-higgsino-higgsino interaction 
terms also produce contributions to chargino and neutralino mass matrices and mixings. 

Completing the program of Antoniadis et al [12] classified higher-dimensional 
operators including the effects of general field redefinitions and loop corrections. We will 
not include any of their operators beyond those of [7] ; this will be explained further below. 

Shortly before this work was concluded, Cheung et al considered effective operators 
in the light Higgsino scenario [13]. We will comment on the relation between their work 
and our work below. 

2 Review of the MSSM Lagrangian 

This section is a very brief review of the MSSM Lagrangian, in order to set notation. Our 
conventions are summarized in appendix |Xj 

2.1 The MSSM superpotential 

The MSSM superpotential is given by: 

VTmssm = -uyuQHu + dydQHd + eyeLHd - fiHuHa. (1) 

where Q are the left-handed quark superfields, u and d are right-handed quark super- 
fields, Hu and Hd are the Higgs superfields, L are left-handed lepton superfields, e are 
right-handed lepton superfields, and yM,yd,ye are Yukawa coupling matrices. We have 
suppressed family indices and SU{2) contractions. (We use the conventions of Haber and 
Kane [13], that are summarized in appendix [A}) 
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The soft supersymmetry breaking Lagrangian is (except for the soft terms in the Higgs 
potential, that we write in Q below): 

Vsoft = {-eAeYeLRd- dAdydQHd + uAuyuQHu + h..c.) 

+ Q^m^qQ + Dm\L + u^mlu + d)mld + e"^m^e (2) 
+ ^ [MiBB + M2{W^W'^ + 2W+W~) + Ms^^ + h.c.) 

where the A are general soft trilinear terms, the m are general soft sfermion masses. 
Ml, M2, M3 are gaugino masses, and the tilded fields are superpartners of Standard Model 
fields. In the DarkSUSY introduction paper p^, the SU{2) contraction is the opposite of 
Haber and Kane, so we have converted eq. (|2| to our conventions (again, see appendix [A| 
for some more details). 



2.2 The MSSM-7 

For later phenomenological purposes we will truncate the 124 MSSM parameters (the 
"MSSM-124") to seven parameters, the MSSM-7. These seven parameters are given at 
the weak scale. They are: a single gaugino mass scale M2, a single universal scalar 
mass mo for squarks and sleptons (replacing the general m matrices in ([2])), the two 
trilinear soft terms At and Ai, (replacing the general A matrices in ([2]) by Ag = 0, A^^ = 
diag(0, 0, Afo), A„ = diag(0, 0, A^)), the ratio of Higgs vacuum expectation values tan/3, 
the supersymmetric Higgs mass yU and the CP-odd Higgs scalar mass itiao- The last three 
parameters will be discussed in more detail in the next section. 

To reduce the three gaugino masses Mi, M2 and M3 to a single scale M2, we impose 
the gaugino mass unification condition 

M3 = — sin^ ewM2 = - — cos^ ^wMi (3) 
a a 

so we can choose M2 as the free parameter. From the effective field theory perspective 
we adopt in this paper, we have no real reason to expect standard gauge unification to 
hold, but since we want to compare to existing models, we will simply focus on this line 
through (Mi,M2,M3) space for simplicity. 

2.3 The MSSM Higgs potential 

The tree-level MSSM Higgs potential is given by 

Vmssm = + mjjjHlH^ + + ml jH^.Hd - {bH^Hd + h.c.) 

(HlH^ + HlHdf - A{HMHHuHd)] + ^{HIH^ - HlHdf , (4) 
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and after writing out SU (2) contractions (using the conventions in appendix |A]) and 
rearranging terms: 

+ {bWH^ - KH'd) + h.c.) 

+ \{9' + 9'){\Hy + \Hi?-\H'A'-\H-\^Y + \g^HtH^^^ , (5) 



which reproduces the expression in [16]. It will be relevant in the following that the 
MSSM quartic Higgs terms are exactly supersymmetric (i.e. completely determined by 
the couplings g and g' of the ordinary electroweak theory, and they are in turn determined 
by experiment). The reason that the quartic Higgs terms are exactly supersymmetric is 
simply by construction: the MSSM only includes soft supersymmetry breaking operators, 
which by definition are of dimension three or lower. In particular, the quartic Higgs 
coefficients (let us schematically call them A) are quite small, and any additional small 
contributions from effective operators may have important effects. For example, in eq. 
Q we have at the Z mass 

A ~ ^ J ^ 0.07 . (6) 
8 

The parameters of the Higgs potential are restricted such that neutral components of 
the Higgs doublets and receive real and positive vaccum expectation values 

Vu={H:), Vd={H'd)- (7) 

At tree level, the vacuum expectation values are related to the and boson masses 
by 

v' = vl + vl = = ^ - (174 GeV)^ (8) 

g2 g,2 g2 

thus only the ratio between Vu and Vd is unspecified. This ratio is denoted 

tan/?=— . (9) 

Vd 

In terms of the Lagrangian parameters {mH^,mH^, fi, b} and the couphng constants g, g', 
the vacuum expectation values are determined by 

2b 

sin 2/3 = . — — (10) 



9^ + 9" Wl - sin^ 2/5 '''' ' ^ ' 
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The expectation value v is fixed to the experimental value eq. ([8|, so the point of equations 



(10), (11) is simply that the original 4 parameters {mH^^mH^, fi,b} in the potential in 



eq. ([5]) can be traded for 3 parameters {tan /3, /x, b} using one constraint We will also 
trade b for m^o below. 

The Higgs fields are expanded around their vacuum expectation values 



and 



Rf-. 



G+ 



We diagonalize the quadratic part of the Higgs potential by the rotation matrices 



Rn 



cos a sm a 



Rpo 



sin cos 
-cos/5° sin/30^' 



Rr- 



— sma cos a 
This gives us the following tree-level mass spectrum 



rriQO = = 0, 



sin cos 
— cos sin 13^ 



(12) 



(13) 



(14) 



m 



AO 



(15) 



1)2 + Am\,ml sin^ (2/?) 



AO 



(16) 



where clearly the lighter Higgs has the upper sign, and the heavier Higgs has 
the lower sign. We have now traded the original 4 parameters {mH^jfnH^jb, fi} in the 
potential in eq. ([S]) for the 3 parameters {tan/?, yU, m^o} using the constraint ([S]) and the 



trading equations" (10), (11), (15). The charged Higgs fields have masses that are simply 



2 2,2 



(17) 



while for the angles in the rotation matrices (14) we find = = /3, and the angle a 
is determined by 



AO ~t~ Z / /0\ 

sin(2a) = ^ 2~sin(2p), cos(2a) 



2 2~ 



COS (2/3). 



Now that we have the Higgs sector masses and mixings, we can read off from eq. (16) 



the familiar statement that the tree-level lightest Higgs mass satisfies rriho < rriz- This 
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is the Higgs version of the "httle hierarchy problem", that fairly large loop corrections 
are needed to bring the Higgs mass above the LEP bound 171^0 > 114 GeV. (See section 
(§ for some more details on how we implement this experimental bound, including loop 
corrections.) The little hierarchy problem has been studied for a long time (see e.g. [17j) 
and was the motivation for [71 El [9] to consider the corrections we are about to include. 

Why is the tree-level Higgs mass so low in the MSSM? Let us consider the toy 
potential 

Vtoy = + (19) 

which has a minimum at </)^j„ = f ^ = yU^/(2A), and the mass squared at the minimum is 
= V"{v) = Ajj?, that is 



m = 2V2Xv 

< 2V2 ■ 0.07 ■ 174 GeV ^ 120 GeV , (20) 

using the rough value of the quartic coupling A from (|6|. As emphasized above, the MSSM 
quartic Higgs coupling A is small due to the exact supersymmetry of MSSM dimension 
four operators. It is this smallness in eq. (|6| that we want to relax a little, so that the 



mass in eq. (20) can be somewhat larger. 



2.4 MSSM neutralinos and charginos 

The MSSM spectrum also includes neutralinos and charginos. We will not repeat the 



standard discussion here, but only give results in the modified theory in sections |4.4| and 
4.5 , except for two things: first, let us already note that we impose that the neutralino Xi is 



the Lightest Supersymmetric Particle (LSP), so when we write ttilsp it always means m^o. 
The mass eigenstate is found from rotating gauge eigenstates ip'^ = (5°, W^, H^, H^Y^ as 

X? = iVii5° + iVi2l^o + iVi3iJ,° + iVi4^° . (21) 

Second, we define the gaugino fraction of the LSP in terms of the matrix elements A'^ij as 

Zg= |iVnP + |iVi2r (22) 



As is clear from (21), this simply says how much of the LSP is gaugino (i.e. or W^). 

When presenting numerical results, we typically prefer to use the physical parameters 
(Zg,mLgp) rather than model parameters like (/i,M2). To provide an idea of the relation 
between the gaugino fraction and the model parameters, it is useful to note that (see e.g. 

m) 

I - Zg sin^ Ow f — 1 when < Mi < /i (23) 
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and 

\ 2 



1 ± sin 2/3 / sin^ cos^ 



Zg — - " + —-pf- m'z when < /i < Mi . (24) 



Here the sign is minus for the neutrahno we call A^3 in eq. (62) below, and plus for NA. 
We impose the gaugino mass unification condition ([s]), so we can always express the Mi 
in these equations in terms of M2. Thus if we set 

mz < M2 < /i = 1000 GeV (25) 

we find 1 — Zg ~ 10^^ (mostly gaugino), and if we set 

< /i < M2 = 1000 GeV (26) 

we obtain ~ 10~^ (mostly higgsino). This is somewhat representative of certain 
parameter sets we will use later, e.g. table [TJ 

This concludes our short review of the MSSM. For more details, we refer to review 
articles such as 1161. 



3 Effective theory 

We now consider allowing for higher-dimensional operators in the Lagrangian of the 
MSSM, with small free parameters as coefficients. For more on how those parameters 
can be related to parameters of specific underlying theories, see appendix [Bj To be pre- 
cise, we will incorporate operators of effective dimension > 4 (counting the power of 
1/M only), but scaling (naive) dimension > 4 (counting the total mass dimension of the 
coefficient, including M but also /i and the scale of supersymmetry breaking msusv)- 

To be useful, the effective theory should not add too many new parameters. The 
number of effective operators one could a priori consider adding to the MSSM at scaling 
dimension four and five is in the several hundreds. We will restrict consideration to 

• operators that preserve baryon and lepton number, which must be a good approxi- 
mation. 

• dropping the fairly large number of operators that can be added in the squark sector 
|19j|^ This may or may not be a good approximation, but is only a simplifying 
assumption. To seriously study physics of effective operators in the squark sector, 
we should consider more general squark couplings and masses than the MSSM- 7 we 
use in the parameter scan below, which has universal masses and couplings at the 
electroweak scale. 

■^These include certain operators that have recently been relevant in a different dark; matter context 
in [50] , based on general gauge mediation as discussed in [5TJ [55] . 
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• dropping CP-violating higher-dimensional operators. To go beyond this, we should 
again go to a more general starting point than the MSSM-7 we use in the parameter 
scan below, which already has almost all the potential CP violation in the MSSM 
set to zero in order to easily pass experimental constraints. 

With this combination of approximations and simplifications, it was pointed out by 
that there are only two dimension-five operators that can be added in the superpotential 
in the Higgs sector, and a reasonably small number of dimension-six operators in the 
Kahler potential. In this paper, we consider the mass scale of new physics M to be 
sufficiently large that the dimension-six (1/M^) operators do not contribute appreciably, 
e.g. M ~ 5 - 10 TeV. 

Consider the BMSSM superpotential 

W = H^MSSM + (27) 

where Wmssm is given in eq. ([T|, and is the following operator at effective dimension 
fiveS 



M 



{2i 



Here A is a dimensionless number and M is a large mass parameter, though M -C Mqut- 
Supersymmetry breaking operators can be parametrized by introducing a coupling in (28) 
that depends on a spurion field Z as 



A^A(1 + Z) 



(29) 



where msusY is the scale of supersymmetry breaking. 

The new terms in the Lagrangian are now calculated (see appendix lAl) from the su- 



perpotential W in eq. (27). From the supersymmetry-preserving part (28) we find the 
following contribution to the Higgs scalar potential, for any scalar field (pf. 



i 

E 



dW 



MSSM 



d(j)i d(f)i 



+ 



MSSM 



+ h.c.+ C(l/M^ 



(30) 



The cross term gives the 1/M operator we are interested in. Explicitly, from the super- 
symmetry preserving operator in (28) we obtain the following contribution to the Higgs 
scalar potential 



5^1 = 2er{H^Hd){HlHu + i/Ji^d) + h.c. 



(31) 



^This is not the same A as in the previous section; this A only occurs in this section and in the 
Appendix. 
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where we introduced the dimensionless number 

^1 = ~^-- (32) 

Notice that ei is dimensionless and involves both fi and M, since it comes from the cross 
term in (30). Thus, (31) has effective dimension five and (obviously) scaling dimension 
four. 



Similarly, from the supersymmetry breaking operator in (29) we get the following 
"soft" contribution to the Higgs scalar potential 

SV2 = e2iHuHd)^ + h.c. . (33) 

Here we introduced the dimensionless number 

£2 = -^msusY- (34) 

Although we called the contribution 6V2 "soft" , it is clearly not soft in the sense of being 
dimension three or less. Martin [TT] calls it technically hard. The point of the classification 
into hard and soft in the MSSM is that hard terms could spoil the relatively nice UV 
behavior of the MSSM. However, since 6V2 is only an effective theory contribution, that 
will be matched to an underlying theory at the scale M (see appendix [B]), it will not affect 
the UV behavior in any important way. 

It is already clear at this point that when we give expectation values to the Higgs 
fields, the operators 6Vi and 6V2 will affect the minimization of the Higgs potential and 
could potentially affect the mass of the lightest Higgs field. 

Finally, using (83) in appendix |A| the supersymmetry-preserving operator (28) also 
produces the following Higgs-Higgs-higgsino-higgsino interaction^ 

= +^ {2{H^Hd)iHuHd) + 2{HuHa){HuHd) + {H^H^f + (/f„^d)') + h.c. (35) 

These terms introduce new interactions, and when the Higgs fields acquire vacuum ex- 
pectation values, also new mass terms for neutralinos and charginos. 

In general, there will also be dimension-five cross terms of (28) with Yukawa couplings 
in the MSSM superpotential ([T]), producing for example 

l-.. + Qg + ^^T (36) 

where H is some trilinear combination of Higgs fields. However, these only give contribu- 
tions to the additional dimension-five operators in the squark sector, that we have already 
set to zero. Further, the higher-dimension terms do not contribute directly to Yukawa 
couplings. They can contribute to processes with 3-particle final states, but we neglect 
those processes. 



^Here we disagree with [7, eq. 18] on the overall sign, since our ei is defined by eq. (32 1. See also 
appendix |Xj 
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4 Effective theory: final result 



In the previous section we constructed our BMSSM theory. To be precise, we constructed 
an "MSSM5" theory, which is the subset of BMSSM to order 1/M. In addition, we do not 
consider the most general MSSM5 theory, but imposed further restrictions on the 1/M 
operators, as detailed in the beginning of the previous section. With this understanding. 



we will consider the MSSM Lagrangian ^ plus the terms (31), (33), (35) 



= -6Vi-6V2 + 6^3 

= - (2e^{HMiHlH^ + H^M + h.c.) - [e^iHM' + h.c.) (37) 



+ 



/i* 



(2iHMiHuH,) + 2iHMiHM + {H^H.f + {H^H^f'^ + h.c. 



As before, our conventions are spelled out in appendix |Aj The terms in the second 
line change the Higgs potential, and those in the third line add new Higgs-Higgs-higgsino- 
higgsino couplings at scaling dimension five, and affect the neutralino and chargino masses, 
mixings and couplings to Higgs bosons. Thus, we have two new free model parameters: the 
small dimensionless numbers ei and €2- Our restriction to CP-conserving new interactions 
(see the previous section) lets us restrict ei and €2 to be real, but they can be positive or 
negative. See also section |8] for some comments on more general BMSSM models. 

As long as we stay at energies sufficiently below the energy scale M, it is a standard 
exercise to calculate cross sections including also the nonrenormalizable interactions in 



the second line of (37). For a brief review of how this is done, see appendix [B| The 
Feynman rules due to the new operators are presented in appendix [C| 

We now proceed to see how the corrections affect the Higgs potential, the neutralino 
masses and mixings, and the chargino masses and mixings. These different aspects of the 
BMSSM corrections can all affect the relic density of dark matter in different regions of 
parameter space, as is investigated in section [7| 



4.1 BMSSM Higgs potential 



Now we turn on ei and €2 from (37) to small but nonzero values. The minimum expressed 



in terms of the vacuum expectation values f „, Vd will change, so the ratio tan j3 will change. 
However, since the vacuum expectation value v = \/vl^rv^ ~ 174 GeV is fixed to the 
experimental value, we keep v fixed when we introduce the ei 2 corrections. 
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The "trading equations" (10), (11) and (15) are replaced by 
2b 



sin 2/? 



sin 2/3 

2 I 



+ sin 2/5 +2t;^e2 



+ + 8ei sin 2/5 \^ ^1 - sin^ 2/5 

,2 I 2 , ol..|2 



2 2 ol |2 



+ '^Hd + 2|/ir + ei sin 2/? - 2t;^e2 • 



(38) 

(39) 
(40) 



In the effective theory, we use these new equations to ehminate {mu^^mu^^h] in favor 
of {tan/5, f , m^o}. As a consequence, the expressions for the charged Higgs mass and 
the CP-even Higgs mass matrix depend explicitly on ei and 62. (We of course still have 
= for the Goldstone bosons.) For the CP-even Higgs bosons we find the 

mass matrix 



M 



m 



AO 



in terms of sp = sin/?, cp = cos/5, and in the basis {H^^,H'^^) for the real parts of the 
fields. This yields, to first order in ei and €2: 



2 I 2 

m^o + rriz =F 



m 



AO 



- m|)2 + 4m^om| sin^ (2/5) j + <5m^o,^^o (42) 



where 



hO 



HO 



2v^ 62 + 2ei sin (2/5) 



2ei(m^o + "^D sin (2/?) - e2(m^o - m|) cos^ (2/5) 



- m|)2 + Arn\om\ sin^ (2/5) 



2^2 €2 + 2ei sin (2/5) 



2ei(m^o +m|)sin(2/5) -eai 



^0 



m|) cos^ (2/5) 



(m^o — m|)2 + Am\om\?,iv? (2/5) 
The charged Higgs fields are now given masses 

m^± = m\o + + 2t;^e2 . 



(43) 



(44) 



(45) 



For an illustration of eq. (43), see fig. [T] below. 

Since we want to compare BMSSM models with the corresponding MSSM models for 
Ci = £2 = 0, we have to make a choice how to treat the parameters tan /5 and m^o. We will 
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treat them the same way as we treated v, i.e. we assign them the same value as we would 
have in the corresponding ei = €2 = MSSM model. One reason this is useful is the 
following. We know that the partial width of — * qq, which is in principle measurable, 
is given directly in terms of tan j3. Thus, keeping it fixed will prove convenient when we 
study the resonance region (see e.g. fig. [s] below), because the width will not change 
appreciably when we turn on ei_2- In other words, the BMSSM model is physically very 
similar to the corresponding MSSM model in this particular respect, which would not 
have been the case had we let tan/3 vary. Similarly, the pole mass of m^o is a physical 



parameter that we would like to keep fixed. The shift in (40) can then be viewed as a 



finite renormalization, that we absorb in the definition of the pole mass, which therefore 



stays the same. The value for rriAo we input in (42) is therefore the same before as after 
turning on ei 2, and as a check the A^ resonance does not move (see fig. [s] below). 

For the rotation angles we find f3^ = (3^ = f] even with 61^2 7^ 0. That is, the 
functional form of and as functions of /? does not change, although /? changes 
(which we reabsorb in the free parameter tan (3, as noted above). On the other hand, the 
angle a receives a correction, which to first order in ei^2 can be expressed as 

. ,^ Av' cos^(2/?) (2ei(m^o - m|)^ + 62(171% - m%) sin 2/3) 
d)(sm(2a)) = ^ (46) 

{{m\o - m%y + Am\Qm% sin^ (2/?)) 
or, which is sometimes more convenient. 



Av^ sin(2/5) cos2(2/3) (2ei(m^o - m|)2 + e2{m% - m%) sin2/3) 
()(cos(2a)) = ^ ^ . (47) 

{{m\o - m?zf + 4m^o"^| sin^ (2/3)) 



4.2 Effects on the Higgs mass 

One important question is the size of 5mhQ since as previously mentioned, this affects the 
Higgs little hierarchy problem. In figure [T] we show an example of how m^o varies with 
ei and £2, and also how the Higgs mass in the ordinary MSSM depends on loop order. 
The calculation is performed using FeynHiggs [23], EU [25l |26] (see section [s] for further 
details). 

To see why a small quartic coupling A makes the mass shift 5m/jO large even for fairly 



small ei,2, let us consider the toy potential Vtoy of (19) again. If we shift this to 

Koy, shift = -/i V' + A' = A + e (4J 



we find v"^ = yU^/(2A') (compare this to eq. (39)) and = V"{v) = Afi^. Since we want 
to keep V fixed by experimental input, the toy mass = V"{v) will shift: 

Sm'^ = 8ev^ . (49) 
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Hence the relative change in m is 

Thus, for the estimate A ~ 0.07 from ([6]), the relative shift in the toy mass m is of order 
50% for e ~ 0.07. (For the actual BMSSM Higgs mass shifts above, the corresponding 
value is slightly higher than 0.07. )|^ 

In [12], it was argued that to avoid "perturbations" as big as 50%, we should in 
principle impose e ^ 0.07. However, how strictly this should be imposed depends on the 
point of view. We argued that the Higgs mass is unnaturally small due to the smallness of 
A in the MSSM, so the Higgs mass could receive relatively large shifts as we "naturalize" 
the unnaturally small coupling A by small additions in ei,2- In practice, we will allow for 
shifts up to order 20%-30% as displayed in figure [l] This restricts ei to < 0.05. 




-3 -2-10 1 2 3 -3 -2-10 1 2 3 

A,lma AJriia 



Figure 1: Higgs mass. The solid curves all include two-loop corrections. Notice the 
MSSM lightest Higgs mass increases with 1-loop corrections, but decreases somewhat 
when 2-loop corrections are taken into account. For this figure we used M2 = 500 
GeV, /i = 700 GeV, tan/5 = 3, m^o = 500 GeV, tuq = 700 GeV and = At. The 
gray area is the rough estimate m/jO < 114 GeV for the region ruled out by the LEP 
Higgs bound, though we remind the reader that we do not use this bound literally 
(see section l5| . 



4.3 Bounds on e\ and 62 

As mentioned above, one could argue that ei^2 too large will cause unreasonably large 
mass shifts in m^o and hence be inconsistent. We will restrict attention to 171^0 mass 
shifts of at most 20%-30%, which gives 

|ei|<0.05. (51) 



^Pedantically speaking, the actual expansion in (50 1 of course breaks down when e ^ A. We just want 
to illustrate that there are large effects when e ~ A. 
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The 62 correction does affect the Higgs mass less than ei, so in principle one could consider 



turning on a larger value than (51) for €2 as far as the Higgs mass shift is concerned. 



However, as it turns out 62 does not impact our calculations very much, so we will simply 



turn on €2 conservatively as in (51). 



Next in the list of restrictions on ei,2, we have the following estimate [27] on precision 
electroweak bounds on the S and T variables [28] expressed as a bound on the scale of 
new physics M: 



7 



M > 8 TeV . (52) 

If we would apply this bound to estimate allowed values of ei, then for a perturbative 
microscopic theory (see appendix [B| the region of low /i would only support very small 



values ei < 0.01 (see eq. (32)). However, the tension with the bounds on the precision 
electroweak observables S and T comes from effective dimension six (i.e. 1/M^) terms in 
the Kahler potential that provide new interactions between Higgs fields and gauge fields. 
These new interactions come with coefficients called in [7] that are independent of ei^2 in 
the effective theory, and we have set = 0, as stated in section |3| Since our software tools 
do perform checks on precision electroweak observables (see section Isl), we will assume 



that (52) does not need to be strictly applied to our BMSSM models. See also section |8| 
When modifying the Higgs potential, one should check vacuum stability. A priori one 
might think that small corrections can never affect stability. However, the usual derivation 
of the condition that the potential be bounded from below is in a D-flat direction, and 
the MSSM quartic Higgs term is zero in this direction (it is a D-term), so in principle 
ei,2 could affect boundedness from below. Nevertheless, for ei,2 sufficiently small and ei 



nonzero, the potential is still bounded from below since (31) is supersymmetric (Vp in eq. 



(30) is an absolute value squared). The issue then becomes whether to see this, we need 



to keep the effective dimension six operator in (30). For the purposes of this paper, we 
will not investigate this further]^ 



For models with large Higgs mass shifts, i.e. large ei_2 beyond (53) below, there is then 
the related issue of whether the €1^2 corrections could deform the potential so severely 
that a new deeper global minimum develops, and if so how long it will take it to tunnel 
to this new true vacuum. Vaccuum stability under such conditions was checked recently 
by Blum, Delaunay and Hochberg (BDH) in ||22|, who found a criterion to exclude such 
transitions in the BMSSM even for large Higgs mass shifts. Since for our purposes we do 
not particularly focus on obtaining large Higgs mass shifts, we will not impose the BDH 
criterion, though again it would be interesting to study this issue in more detail. 



■^Note that if we prefer, we can always ensure that the quartic term is positive by turning on a 
sufficiently large positive £2, although we did not do so explicitly. We believe our results will remain 
mostly unchanged if we did. 
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To summarize, we will require the following reasonably conservative bounds on the 
BMSSM model parameters ei and 62- 



- 0.05 < ei < 0.05 , -0.05 < £2 < 0.05 



(53) 



4.4 BMSSM neutralino masses and mixings 

In the basis ■0° = {B, W, H^, H^Y , the neutrahno mass part of the Lagrangian is given 

by 



(54) 



where in terms of — sin/3, = cos/? and the usual Standard Model parameters m^, 
sw — sin s-nd — cos Q-w we have 



(55) 



where at tree-level 



f Ml 

M2 



-mzswcp mzswsp \ 
mzcwcp -mzcwsp 
-11 



\ mzswsp 



-mzCwS/3 







(56) 



/ 



(i.e. the same functional form as in the MSSM) and the new couplings introduced above 
introduce the following corrections to the neutralino mass matrix 



6M^o = 



2ei 
II* 



/O 










\ 





2ei 2 

V 

II* 



/O \ 



sin^/? sin 2/5 

\p sin 2/3 cos^/jy 



(57) 



\0 2vuVd vl J 
The mass matrix is diagonalized using a unitary matrix A^, which is such that 

N*A4^oN^ — diag(m^o, m^, m^, m^o) 



(58) 



is diagonal and ordered according to m^o < m^o < m^o < m^o. The mass eigenstates are 
given by 



(59) 



In the limit where the electroweak symmetry breaking terms can be considered to be 
small, i.e. 11 or Mi, M2 ^ mz, and for real parameters, the eigenvalues of the neutralino 
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mass matrix are given by (see e.g. [18 j ) 



M| sin^ OwiMi + fi sin 2f3) 
Ml - /i2 



mTvi = Ml H _ 2 ^ • • • (^^^^ 



M|cos^gH/(M2 + /isin2/3) 
m^, = M, + + ... (61) 

II ei 2n o • om • ( ^"^l(l + sin2/5)(MlCos2^vy + M2sin2^H/-/i) , 

"^"^ = - - - "S"^^) 2(Mi-;z)(M2-/.) + 

(62) 

I I ^ ei 2n ^ o • om ^ • ^ , m|(l - sin 2(3) (Mi cos^ + sin^ + /^) ^ 

"^"^ = + + '""'^^ + 2(Mi + ;.)(M2 + /.) + 

(63) 

where we have introduced the non-mass-ordered neutrahnos A^l through NA. Note that 
even though the ordering can change due to the ei corrections, this is not a problem for 
the consistency of the theory; the mass corrections due to ei themselves remain reasonably 
small for reasonably large /x, as we show in an example in figure |2] below. 

We include the most significant MSSM loop corrections to the neutralino mass ma- 
trix [291 EHl EH |32], as implemented in DarkSUSY (see section [s]). Although these loop 
corrections are small (on the order of a few GeV), the corrections in the (3,3) and (4,4) 
elements of the neutralino mass matrix can be important, since those elements are zero 
in the MSSM at tree-level. The BMSSM corrections also contribute to those elements of 



the neutralino mass matrix (see eq. (57)), so it is conceptually important to include loop 
corrections in the BMSSM as well, since otherwise the effect of ei,2 would appear greater 
than it really is. The loop corrections are themselves in principle affected by the BMSSM 
operators, but this effect is higher order and is neglected here. 



4.5 BMSSM chargino masses and mixings 

In the basis ip± = {ip^ ,ip~) where ip^ = {W~^, -f^J')'^, ip~ = iW~ , H^)'^, the chargino mass 
terms in the Lagrangian are given by 

^3-^(^±)^^x+^±+ h.c. 

= -^{^+fXi^. - ^{^.fX^^+ + h.c. (64) 

where the mass matrix A^3f+ takes the following block-diagonal form 

M,. ^ (» f ) (65) 
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with 



where at tree level 
X 



M2 gvu 
gvd 



X = X + 6X 



M2 V2Mw sin 

\/2Mw cos l3 /i 



(66) 



(67) 



(i.e. the same functional form as in the MSSM) and the dimension-five operators intro- 
duced in eq. (37) give the following correction to the chargino mass matrix: 

6X 



— f ^ sin 2/3 
/i* 




1 



(6^ 



The matrices X and X'^ are diagonalized using two unitary matrices U and V such that 
U*XV^ = diag(m^„ m^J, V*X''U^ = diag{m^,,m^,) (69) 
with the chargino masses on the diagonal, and where the mass eigenstates are given by 



For the standard MSSM case with ei 2 = one finds 



2 

w 



(70) 
(71) 



T ^J{\M2\'' + W + 2m2,)2 - 4|^M2 - ml, sin2/5|2) 



(72) 



For future reference, we define x (without the index 1, 2) to be the lightest chargino. Now, 
we see from eqs. (67) and (68) that the corrected chargino masses are simply obtained by 

(73) 



replacing /x in (72) with 



H H — ^-v^ sin 2/3 
11* 



(Note that we do not use this /i' anywhere else.) In the case of small real we have 



m~± = I u| ( 1 H — sin 2/? -|- . . . 
^2 \ 



(74) 



In particular, we see that the chargino-neutralino mass splitting receives an ei-dependent 
correction. We plot an example in fig. |2} In the higgsino region, we can write a simple 
expression for the neutralino mass crossing (which incidentally does not need to be 
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Figure 2: Masses and mass crossing for the StHelena^^-* model (see tables [l] and [i] 
below) for e2 = and varying ei. The rough estimate mj^o > 114 GeV of the LEP 
Higgs bound is shown, though we remind the reader that we do not use this bound 
literally (see section p|. The point ei = e'i°^^ marks the crossing point, which is 



given in eq. (75 ). With the rough estimate of the LEP Higgs bound, the point e'^'^^^ 
is in the excluded region in this example. For further discussion, see the main text. 



at the same place as the neutralino-chargino mass crossing, despite appearances for our 
benchmark model in fig. |2|): 



- /xsin2/?(s^M2 + cl^Ml - /i^) + ^?{sl,Mi + c^Ma) - MiM2(s^M2 + c^Mi)) . 
Then the neutralino-chargino mass splitting is a simple linear function of ei: 

I ^x± ~ ^NA II ei < Cl 
where m tvs and m ^4 are the expressions in ( 62 ) and ([63]) 



*^In 13J, it appears that when ei increases, the Higgs mass increases while the chargino mass decreases, 
the opposite of figure[2] This is due to the use of two incompatible conventions for ei in their expressions. 
Since this probably affects the results for the relic density in that paper, we find it difficult to compare 



their results to ours. See also section 7.3 
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5 Computational tools 



For every BMSSM model parameter set we consider, we calculate the resulting particle 
spectrum and relic density using DarkSUSY [T3], a publicly-available comprehensive nu- 
merical package for neutralino dark matter calculations. 

Starting from the DarkSUSY-5 .0.4 version, we have extended the code to also include 
the BMSSM corrections due to the effective terms (37). The neutralino and chargino 
mixing matrices of eqs. (57) and (68), as well as the three-point vertex corrections and 
the new four-point vertices of appendix [C] have all been fully implemented into the codej^ 

Within the DarkSUSY package the Higgs boson spectrum is determined at two-loop 
level using the external FeynHiggs-2 .6.4 program [221 [211 [2S1 12S]- We have implemented 
the Higgs boson mass corrections of eqs. (43), (44) and (45), the rotation angles (46), 
(47), and the neutralino and chargino mixing matrix corrections in eqs. (57) and (68) in 
FeynHiggs-2 . 6 .4. We do not include additional radiative corrections due to the three- 
point vertex corrections or the new four-point couphngs, but we expect our results to be 
largely insensitive to this deficiency. In addition to the masses, we also extract the total 
and partial widths of the Higgs bosons from FeynHiggs. 

It is also important in this context to check the LEP accelerator bounds on the Higgs 
boson masses carefully. In the BMSSM (or any other theory with a Higgs sector different 
from the Standard Model Higgs boson — including the MSSM) the Standard Model lower 
Higgs mass limit of 114.4 GeV strictly speaking does not apply; in principle, the bound 
should be determined model by model from experimental data. One could worry that if 
the bound would effectively be weakened when going from the MSSM to the BMSSM, 
that would take away some of the motivation for including the corrections (43) in the first 
place. 

For this reason we call HiggsBounds-1 . . 3 from DarkSUSY. The HiggsBounds 

code takes actual Higgs boson production cross sections and partial widths as input, 
computes rates into the important LEP search channels, and determines whether or not 
the model is excluded by comparing the rates with existing LEP data. From our parameter 
scans we note that the simple LEP bound rriho ^114 GeV in general is respected to good 
accuracy within the BMSSM (including the MSSM), and the simple interpretation of the 
BMSSM corrections as providing new models by raising m/jO works wellj^ 

Finally, we note that DarkSUSY imposes many additional accelerator constraints on 
the models, including bounds on sparticle masses, the rate of b — >• S7 and the electroweak 
observable p. 



^We intend to make this extension of DarkSUSY public in the future, though perhaps not in the near 
future. 

^"It is useful to keep in mind that for special parameter sets, the LEP bound 114.4 GeV can be lowered 
significantly |28j — this is a manifestation of the model-dependence of the bound we emphasized in the 



paragraph above. 
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6 Relic density: strategy 



We explore the features of the BMSSM by scanning extensively over the full nine-dimen- 
sional parameter space, i.e. over the parameters M2, /i, tan/3, m^o, mo. At, A^, ei and €2- 
(For a reminder of what these parameters are, see section |2}) The scans are carried out 
with the help of DarkSUSY [T3], in which we have implemented the BMSSM as described 
in section |5] above. 

For each model we calculate quantities like mass spectra and relic density, including 
coannihilations (see e.g. [S])- We check them against various accelerator constraints, 
including bounds on sparticle masses, Higgs boson masses, the rate of b — > S7 and the 
electroweak observable p. 

Of particular importance for the BMSSM are limits on the lightest Higgs boson and 
chargino. To determine the former we use HiggsBounds [33] (see section [s]), while for the 
chargino bound we adopt 

m^± > 94 GeV, (77) 



which is the current standard lower mass limit [28]. (Recall from section 4.5 that we define 
X"^ without any 1, 2 subscript as the lightest chargino.) In close analogy to the discussion 



of the Higgs mass bound in section [5j above, the chargino bound (77) is somewhat model 
dependent, and should not be thought of as applying with high precision for an arbitrary 
model. Depending on particular properties of a given model, the limit may move several 
GeV in either direction [28], and it is not simple to implement an accurate bound even 
on a model by model basisj^ For illustration purposes, we simply use a sharp 94 GeV 
chargino mass bound as in eq. ([77]), with the understanding that this is not very precise. 

Since the BMSSM is only valid below the scale of new physics, which we take to be 
M ~ 5 — 10 TeV (see section [3]), we also impose a maximum cut on sparticle masses by 
requiring IM2I, |mo| < 2 TeV. This cut could be raised if we raise M, but then our 



new parameters ei, 62 will typically be rather small (see eqs. (32), (34)). 

For the Lightest Supersymmetric Particle (LSP) to be a viable dark matter candidate 
we require it to be a neutralino, and that it provides a relic density in agreement with 
the value measured by WMAP [35j. To be precise, we use the constraint 

n^/i^ = 0.1099 ± 2 ■ 0.0062, (78) 

which is the 2a result when combining ACDM, Sunyaev-Zeldovich and lensing datasets 
[3S] . Here fl^ is the dark matter relic density as a fraction of the critical density, and the 
dimensionless parameter h is the Hubble constant in units of 100 km Mpc~^s~^. 



^^By analogy with the HiggsBounds package, it would be most useful if someone wrote a 
CharginoBounds package! 
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In addition to our broad parameter scans, we also highlight regions of particular in- 
terest for the BMSSM by investigating details of a number of focused scans and certain 
benchmark points. 

We will show most of our plots in terms of the physical parameters mLSP and the 
gaugino fraction Zg. See the discussion in section 2.4 for some details on the relation 
between Zg and the model parameters M2 and /i. 



7 Relic density: results 

7.1 General parameter space scan 

We have scanned the MSSM parameter space between generous bounds allowing mass 
parameters up to several TeV, At and Ab between — 3mo and 3mo and tan/3 between 
1 and 60. Our sample of MSSM models is larger and we have about 110 000 models 



that pass all constraints (including the WMAP constraint on il^h of eq. (78)). For the 
BMSSM, we have about 11000 models that pass all constraints. To be clear, the fact 
that we have ten times more MSSM models than BMSSM models does not mean that it 
was in any sense harder to find BMSSM models than MSSM models (of course for 61^2 
infinitesimally small, the models tend to be indistinguishable), it merely reflects the fact 
that at the time of writing, we have performed more MSSM scans. 

In figure [3] we show the results from this scan projected onto the (mLgp, Zg/(1 — Zg)) 
plane, where mLgp is the mass of the lightest neutralino (our dark matter candidate) and 



Zg is the gaugino fraction of the LSP, as defined in eq. (22). The blue circles and red 
dots in the figure both represent regions in the (mLSP,Zg/(l — Zg)) plane where we find 
models consistent with all imposed accelerator constraints and that provide the correct 



relic abundance of eq. (78). The red dots correspond to regions where we only find models 
with ei 7^ and/or €2 7^ 0, i.e. regions which the ordinary MSSM is unable to reach. Since 
the set of all MSSM models in figure |3] bears a vague resemblance to a map of Africa, we 
label various interesting regions by African countries (and one island) in the appropriate 
locations. 

The only feature of the BMSSM that is immediately obvious in figure |3] is the existence 
of completely new models with light (mj^gp ~ 75 GeV) higgsino-like neutralinos; the red 
dots around the regions marked "St Helena" and "Ghana". We show a plot zoomed in 
on this region in fig. |4| and elaborate on this in sections 7^ and 7^, respectively. 



The BMSSM can also by construction naturally accommodate light top squarks with- 
out tension with the mass bound on the lightest Higgs bosons, unlike the ordinary MSSM 
(see section 4.2 above for more details on this so-called "little hierarchy problem"). This 



property is investigated in the dark matter context for the "Niger" region in section 7.4 



We note that these light-stop BMSSM models, and many of the other ones we study be- 
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Figure 3: MSSM and BMSSM models that pass all accelerator constraints and fur- 
nish the WMAP relic density. Shown in the figure are the locations of our benchmark 



scans of sections 7.3 (St Helena), 7.4 (Niger), 7.5 (Mauritania) and 7.6 (Ghana). 



low, are not visible in the commonly used "projected" parameter plot in fig. [3} i.e. when 
projected onto the (mLSP, Zg/{1 — Zg)) plane they are covered by MSSM models. Thus, 
the fact that fig. [3] is dominated by blue circles should not be taken to mean that new 
BMSSM physics occurs only where the red dots are located. This will be discussed further 
in the following sections. 



In the remaining sections 7.3 through 7.6 we investigate these more subtle aspects of 
BMSSM models through a few selected benchmark scans. 



7.2 A slice through parameter space 

Figure [3] is useful for hints for experimental searches by showing what the genuinely new 
BMSSM models in the (mLSP, Zg/ (1 — Zg)) plane are, but it does not express much of the 



25 



i 1 



10 



Ghana 



Berg, Edsjo, Gondola, Lundstrom and Sjors, 2009 

oooooooooooooooooooo 
ooooooooooooooo 

(OOOOOOOOOOO - 
(•OOOOOOOOO* = 

••••oooooo* 
•••ooooo* 

• • • • o o o • • ,.--,--> 

• •••00«« - 



MSSM + BMSSM 
BMSSM only 



10 



St Helena'' 
St Helena* 



,<-) 



o 



•'• 

o 
o 
• 

o 

« — 



• • • 

• o 



10 



50 



60 70 80 90 100 

Neutralino Mass (GeV) 



Figure 4: Zoom- in of the bottom-left region of figure k3l Shown in the figure are 



our benchmark scans of sections 7.3 (St Helena) and 7.6 (Ghana). In each group 
of curves, the red curves are positive ei, green curves are negative ei, and the blue 
curves are the corresponding MSSM models (ei,2 = 0). 



dark matter physics of individual classes of models. To investigate this in more detail, 
we consider a specific slice of the full BMSSM parameter space. For this slice, we fix all 
parameters but M2 and /i according to tan /? = 15, m^io = 500 GeV, mo = 1000 GeV, 
At = A^ = 500 GeV, 62 = 0. 

The result is shown in fig. [5} where the black and red points correspond to models with 
ei = and e\ = 0.05, respectively, that are consistent with accelerator constraints and 
provide a dark matter relic density in agreement with the WMAP 2a bounds in eq. (78). 
Note that fig. |3]is the projection of many different slices similar to the left panel of fig. [5] 
on top of each other. (The position of the particular slice in fig. |5] in the full scan of fig. 
[3] can be understood from the labelled benchmark scans.) This should be interpreted as 
follows. Although models in the red BMSSM regions of fig. [5] that are not on top of some 
corresponding MSSM models in black may seem to provide physically "new" models, this 
is not always so, since we might be able to reach physically similar points by varying the 
seven parameters in the MSSM-7 (i.e. leaving this slice) instead of varying the BMSSM 
parameters €1^2- 
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Figure 5: Slices through parameter space, in the (mLSP, -^g/(l — Zg)) plane and 
in the (M2,/i) plane. The red points are ei = 0.05, and the black points are ei = 0. 
Here 62 = 0. The blue points are ei = 0.05, but with \fi\ > 2 TeV. Shown in the 

(St Helena(+)), 



7.3 



7.5 



figure are the locations of our benchmark scans of sections 
(Mauritania) and |7.6| (Ghana) . In each group of benchmark scans, the red curves are 
positive ei, green curves are negative ei, and the black curves are the corresponding 
MSSM models (ei,2 = 0). 



But sometimes, there is something that can only be achieved in the MSSM-7 by varying 
a few specific parameters. For example, near the slice through parameter space shown 
in fig. [5} the position and width of the resonance in fig. [s] is fixed for fixed rriAo and 
tan/5 (see section 4.1). In the MSSM slice (black), this means the tree-level Higgs mass 
is then also fixed, and the resonance on the left (in the "Mauritania" region) can only 
move by adjusting the loop corrections through rriQ and At. In the BMSSM, the position 
of the h'^ resonance can be adjusted by adjusting ei 2, as can clearly be seen in the figure 
(the motion of the peak to the right when we turn on ei). So, we do expect physically 
"new" models in the BMSSM case, even though they appear on top of each other in the 
projected figure [sj We investigate the motion of the peak in the Mauritania benchmark 
scan in section 17.51 

The differences between the MSSM and the BMSSM around the mixed gaugino- 



higgsino region labelled "Ghana" are explored in section 7.6 



For convenience we also express the parameter space slice of this section in terms of 
the input parameters M2 and /j,. The result, which is a deformed version of the left panel 
in fig. [5} is shown in the right panel of fig. [5] 
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Parameter 


St Helena(-) 


St Helena(+) 


Ghana^"'"-' 


Mauritania 


Niger 


M2 


1000 


1000 


300 


120.5 


500 




-101 


94 


118 


500 


700 


tan P 


10 


15 


15 


15 


3 


m^o 


500 


500 


500 


500 


1000 


mo 


1000 


1000 


1000 


1000 


298 


At 


2000 


500 


500 


500 





At 


2000 


500 


500 


500 





ei 


-0.05 


0.05 


0.05 


0.05 


0.05 


£2 














0.05 



Table 1: Benchmark models, as referred to in the main text. All parameters with 
mass dimension are given in GeV. 





St Helena(-) 


St Helena(+) 


Ghana(+) 


Mauritania 


Niger 


"^LSP 


78 


76 


79 


60 


248 




8.1 • 10-3 


1.1 • 10-2 


2.5- 10-1 


1.1 • 102 


1.2 ■ 102 


mr-0 
A2 


115 


118 


143 


116 


478 


m^o 


505 


505 


165 


509 


706 




1006 


1007 


325 


510 


727 


m^± 


99 


95 


108 


116 


478 


rrih 


117 


122 


125 


124 


130 


mi. 


825 


1055 


971 


973 


276 




1173 


971 


1055 


1053 


397 




0.114 


0.108 


0.115 


0.105 


0.111 



Table 2: Spectra and relic density for the benchmark models of table [1} All 
parameters with mass dimension are given in GeV. 

7.3 Light higgsino LSP (St Helena) 

To obtain a dark matter relic density in the region favored by WMAP, the light higgsino 
region in the MSSM is problematic for several reasons. The first challenge is that whenever 
"^LSP ^ ~ 80 GeV, annihilation of LSPs into pairs of W bosons is effective and pushes 
the relic density far below the value favored by WMAP, so higgsino-like neutralinos can 
typically provide good dark matter candidates only if they are sufficiently heavy (> 900 
GeV). But what about neutralinos with masses mLSP ^ 80 GeV below the W threshold, 
where the annihilation cross section should be significantly lower? Here we instead run 
into two other problems, as we now describe. 

Firstly, the 94 GeV lower limit on the chargino mass requires the mass splitting Am± = 



m^± — mLSP between the lightest chargino and neutralino (see eq. (76)) to be at least 14 
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StHelena(-) 


St Helena(+) 


Ghana^ •* 


Gliana(+) 


Mauritania 


Niger 


BMSSM models 




M2 


1000 
-86 ...-121 


1000 
79 ... 116 


300 
65 ...218 


300 
88 . . . 226 


102 .. . 132 

500 


485 ...515 
700 


Corresponding MSSM models (ei = €2 = 0) 




M2 


1000 
-64 . . . -104 


1000 
67 ... 106 


300 
72 ... 220 


300 
72 ... 220 


102 .. . 132 

500 


485 ...515 
700 



Table 3: Benchmark scans, as referred to in the main text. Ah parameters with 
mass dimension are given in GeV. 



GeV when mLSP ^ 80 GeV. However, in the higgsino region, the mass scales of both the 
lightest (and next to lightest) neutralino and the lightest chargino are set by the same 
value /X. Typically the particles are nearly degenerate in mass, which means that all such 
models will be excluded by the chargino mass bound. 

Secondly, even if the chargino bound could be relaxed somewhat in particular models 
(for example for Am± < 3 GeV ^28j , see also section [6] for further comments on the 
chargino bound), chargino coannihilation is very important at such small Am± and causes 
the relic density to come out low also when m-LSp < mw 

To summarize, it is generally difficult to find light higgsino LSPs in the MSSM with 
the relic density favored by WMAP. 

This conclusion can be substantially altered in the BMSSM. As was discussed in section 
4.5 and shown in fig. |2| the ei parameter can introduce a larger mass splitting Am±, and 



hence avoid both the chargino mass bound and the coannihilations at the same time. This 
is the reason for the existence of new models with higgsino-like neutralinos just below the 
W threshold in figures [3]j5j 

Note that it is not completely impossible to find valid MSSM models even in this 
region, as can be seen in figures |3] and |4} However, these seem to appear only if parameters 
are finely tuned. The blue circle at the very bottom of the light higgsino-like LSP region 
actually only includes a single MSSM model, which is very special in the following sense. It 
has M2 > fi but not M2 ^ /x, to be precise M2 = 255 GeV, /i = —77.2 GeV, which would 
typically produce a mixed higgsino-gaugino LSP for generic values of the other parameters. 
However, this MSSM model has the specific (small) tan/3 value of tan/3 = 1.16, and then 
the diagonalization of the mass matrix Ai^o (which is given by (56) at tree level in the 
MSSM, though we also include loop corrections) is such that the LSP is higgsino-like 
despite the relatively low value of M2. In the MSSM, mixed or gaugino-like LSPs have no 
problem giving a sufficiently large neutralino-chargino mass splitting, and this particular 
higgsino-like LSP inherits this virtue, so it clears both the WMAP lower bound and the 
chargino mass bound. The BMSSM, on the other hand, naturally provides higgsino-like 
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models with large neutralino-chargino mass splittings without the need for any tuning of 
parameters. This property is apparent in fig. [4], where we see that the BMSSM models 
dominate the higgsino-like region even though the MSSM and BMSSM samples contain 
a similar number of models for the region shown in that figure. 

In table [l] we show the parameters of two benchmark models, St Helena*^"-* and St 
Helena^"''^ with exactly the BMSSM properties discussed above. As can be seen in table 
|2} both of them provide a lightest neutralino mass below the W threshold and relic density 



in agreement with eq. (78), and without violating any accelerator bounds. The models 
differ in the signs of ei and fi, and also come with different values of tan/3. At and Af,. 
The magnitudes of n have been chosen for the models to give the correct relic density. 
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Figure 6: Relic density O^^/i^ along the St Helena^"-' and St Helena^"*"^ benchmark 
scans, respectively (see table [S]). The BMSSM models (ei 7^ 0) are shown in red, 
while the blue curves show the corresponding MSSM models (ei = 0). The WMAP 
2(7 band is shown, as well as the LEP chargino bounds expressed in m-LSP for the 
individual scans. For discussion, see the main text. 



By scanning over values of /i, we produce figure [6| where we plot the relic density 
n^/i^ as a function of the neutralino mass mLSP for the same set of parameters as those 
of the St Helena'-"-' and StHelena*-"^^ benchmark points (except for the value, of course). 
We also show the results for the corresponding MSSM models, i.e. using the same set 
of parameters except for setting ei = (and again scanning over \x). The ranges of \i 
corresponding to the neutralino mass range 60-100 GeV are shown in Table |3} The reach 
of the St Helena(-) and St Helena(+) scans in the (mLSP, ^g/(l ~ ^g)) plane is indicated 
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in figs. [3] and [ij The StHelena*^"^^ scans are part of tlie parameter space slice of fig. [sj and 
lience appear tfiere. 

For eacli individual scan, the m^± > 94 GeV constraint is translated into a lower 
bound on mLSP, as indicated by the jagged vertical lines in fig. [6j As ei is switched 
on, and consequently as Am± is increased, we see how the chargino bound is effectively 
weakened, i.e. moved to the left when expressed in mLSP- The shaded horizontal bands 
display the WMAP favored region of eq. (78). We see that there are models in both the 
St Helena*^"-* and St Helena*-"'"'' scans that can satisfy the chargino and WMAP bounds 
simultaneously. The resulting relic densities for the corresponding MSSM models are 
reduced because of (chargino) coannihilations. We also repeat (see section |6]) that the 
entire parameter range of these scans also satisfy all our remaining imposed accelerator 
constraints, including bounds on the lightest Higgs boson mass and b — > S7. 

The light higgsino LSP region was also studied by |T3] using analytical estimates for 
a few coannihilation processes and using a fixed relative velocity. It is a priori difficult to 
say to what extent this reproduces a full analysis with all relevant processes and actual 
velocity distributions. They bring up a few of the points above, but we find it difficult to 
compare their preliminary analysis to ours. In particular, we cannot reproduce the relic 
density as in their main dark matter result, fig. 4 of [13]. Instead, we find that the models 
with their e > (which corresponds to our ei < to get the same sign as theirs in the 
chargino mass, see footnote [s]) are excluded by bounds on the Higgs mass. 



7.4 Heavy gaugino LSP (Niger) 

Within the MSSM, most models have heavy top squarks ti^2, due to the well-known 
need for large loop corrections to the lightest Higgs boson mass rriho in order for a given 



model to pass current accelerator bounds. (See section 2.3 for a discussion of this.) 
Roughly speaking, needs to be about > 1 TeV to make /i° clear the bound. (See 
section |5] for comments on the meaning of the LEP Higgs bound in models other than the 
nonsupersymmetric Standard Model.) 

Since the new BMSSM parameters ei and €2 introduce corrections to the Higgs boson 



masses according to eq. (43), these can be used to raise the Higgs boson mass even in 
models with light top squarks. This issue has been discussed in the particle phenomenol- 
ogy literature many times (e.g. [H El E]), to which we refer for details beyond those we 
have given in earlier sections. 

An example of a gaugino-like model with light squarks, "Niger", that passes all ac- 
celerator constraints is shown in tables [T] and [2] By adjusting the value of mo, we have 
introduced the right amount of squark coannihilations to find a relic density consistent 



with eq. (78). Note that the corresponding MSSM model, i.e. with the same parameters 



except for ei^2 = 0, is excluded due to its low Higgs boson mass {rrihO ^ 91 GeV). We 
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Figure 7: Relic density Q^/i^ along the Niger benchmark scans (see table [3|. The 
BMSSM models (ei,2 = +0.05) are shown in red, while the blue curve shows the 
corresponding MSSM models (ei,2 = 0). The curves are practically on top of each 
other. The WMAP 2a band is shown, while the LEP chargino bounds appear first 
at much lower masses than shown here. The MSSM models are excluded by the 
LEP Higgs bound for all values of w-lsp in the plot. For discussion, see the main 
text. 



expect this to be a fairly generic phenomenon in BMSSM models, since the alleviated 
"little hierarchy problem" is one of the distinguishing characteristics of the BMSSM. 

In figure [7] we show, by scanning over M2, the relic density as a function of the neu- 
tralino mass, for the same set of parameters as those of the Niger benchmark point (except 
for M2, that is now scanned over). The shaded horizontal band displays the WMAP- 



favored region of eq. ( 78 ) . The figure also shows the resulting relic density for a scan over 
the corresponding MSSM models; the two curves are almost exactly on top of each other. 
This is as expected, simply because all ei corrections in the neutralino/chargino sector 



go as ei//i (see eq. (37)), which is smaller for gaugino-like than for higgsino-like models. 
Except for the Higgs mass bound, which is satisfied for all BMSSM points but none of 
the MSSM points, all models in the scans satisfy the remaining imposed accelerator con- 
straints. The ranges of M2 corresponding to the neutralino mass range 240-255 GeV are 
shown in Table [31 The position of the Niger scans in the (mLSP,^g/(l — Zg)) plane is 



indicated in the overview plot in fig. [3] of section 7A_ 
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Figure 8: Relic density O^/i'^ along the Mauritania benchmark scans (see table 
[3|. The BMSSM models (ei = +0.05) are shown in red, while the blue curve shows 
the corresponding MSSM models (ei = 0). The WMAP 2a band is shown, while the 
LEP chargino bounds appear first at lower masses than shown here. For discussion, 
see the main text. 



7.5 Light gaugino LSP (Mauritania) 

While light gaugino-like neutralinos in general provide relic densities that are too high, 
they could become valid dark matter candidates if the annihilation cross section is in- 
creased by virtue of an s-channel Higgs boson resonance. In fig. [5] we can clearly see both 
the resonance around mLSP ~ 60 GeV and the resonance at mLSP = 500 GeV. Since 
we use the physical m^o as an input parameter in both the MSSM and the BMSSM, the 
position of the A^'-resonance cannot move when we change ei and €2. The resonance 



will move, however, by an amount roughly given by the analytical expression in eq. (43), 
and this is clearly manifest in fig. |5| 

In tables [1] and [2] we give an example of a BMSSM model, "Mauritania", which passes 
all accelerator constraints and for which the resonance brings down the relic density 
into the region favored by WMAP. (Recall that in this region, VtyJi^ tends to come out 
too high and needs to be brought down, whereas in the light higgsino region in section 



7.3l needed to be increased. 



In figure [8] we show the relic density as a function of the neutralino mass rriLSp for 
the same set of parameters as those of the Mauritania benchmark point, apart from M2, 
which is scanned over. We also plot the results for the corresponding MSSM models, i.e. 
with ei = instead, which gives a Higgs boson mass of m/jO ^117 GeV. The ranges of 
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M2 corresponding to the neutralino mass range mLSP = 50-65 GeV are shown in table |3| 
All models in the scans pass our imposed accelerator constraints. The shaded horizontal 
band displays the WMAP bound of eq. (78), as before, and we can see how the neutralino 
masses favored by dark matter constraints slide to the right along with the ei correction 
to the light Higgs boson mass. The asymmetric shape of the curves originates from 
the Boltzmann tail of the velocity distribution in the early Universe. The reach of the 
Mauritania scans in the (mLSP, Z^/il — Z^)) plane is indicated in figures [s] and Isj 



7.6 Mixed higgsino-gaugino LSP (Ghana) 

As can be seen from the separation of the red and black curves in figure |5} the effects 



that give rise to the new higgsino-like models of section |7.3| are present already at quite 
moderate values of Zg/(1 — Zg), i.e. for fairly mixed neutralinos. In tables [l] and |2] we 
give an example of a benchmark model "Ghana^^-*" that provides correct relic density and 
passes all accelerator constraints, in a part of the light mixed neutralino region that the 
MSSM models of figure |5] are not capable of reaching — as long as we stay in the slice 
through parameter space shown in figure |5j See section [8] for some further comments on 
what this means. 

We investigate the mixed region in more detail by scanning over n, for parameter 
values that are otherwise the same as in the Ghana^"*"-* benchmark model. In figure joj 
the relic density is shown first as a function of mLSP (left panel), and then as a function 
of Zg/{1 — Zg) (right panel) in intervals such that the graphs transition smoothly from 
the left into the right panel. We also plot the relic density for the corresponding MSSM 
models (i.e. with ei = 0), as well as for another set of BMSSM models, the Ghana'-"-' 
scan, in which we instead set ei = —0.025. The ranges of fi used in fig. |9] are stated in 
table [sj Except for the chargino bound and (only for high values of Zg/{1 — Zg) in the 
Ghana^") scan) the Higgs bound indicated in the figures, all models satisfy our imposed 
accelerator constraints. The ranges of these scans are indicated in figures |3] and [5j 

The m^± > 94 GeV constraint is translated into a lower bound on rriLSP for each 
individual scan, as is indicated by a jagged line across each separate curve in the left 



panel of figure M (Recall from section 7.3 that as ei is increased, Am± = m^± — mLSP 



is increased, so the bound expressed in /tilsp nioves to the left.) Also shown as a shaded 



horizontal band is the WMAP-favored region of eq. (78). We see that there are models in 



the Ghana'-"''-' scan that can provide correct relic density while at the same time avoiding 



the chargino bound. As pointed out in section 7.3 , this is thanks to the drop in annihilation 



rate right below the W pair production threshold, in combination with larger neutralino- 
chargino mass splitting that makes the coannihilations less important in the BMSSM 
model than in the corresponding MSSM model with ei^2 = (the blue curve). Also note 
how the presence of the m/jO resonance makes the relic density drop sharply for even lower 
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mLSP (GeV) ZJ{\ - Z^) 



Figure 9: Relic density ^x^'^ along the Ghana benchmark scans (see table |3|. The 
BMSSM models (ei 7^ 0) are shown in red and green while the blue curve shows the 
corresponding MSSM models (ei = 0). The WMAP 2a band is shown, as well as the 
LEP chargino bounds expressed in ttt-lsp for the individual scans. The right panel 
begins at Zg = 0.49, which in this slice (see fig. [3]or[5| corresponds to rriLSP = 100 
GeV, where the left panel ends. The Ghana^") benchmark scan (in green) is also 
constrained by the LEP Higgs mass bound, as indicated in the figure. For discussion, 
see the main text. 

(~ 60 GeV) neutralino masses. In the Ghana'-"'"^ scan, however, the chargino bound just 
about rules these models out, but the Higgs resonance is clearly visible in fig. |5]for only 
slightly higher gaugino fractions. Neither of the other two scans shown provide any dark 
matter candidates near the W pair production threshold. The reason that the Ghana*^"-' 
scan (in green) ends abruptly for rriLSP ^ 65 GeV is that the chargino becomes the LSP. 

The right panel of figure |9] shows how the relic density increases, and finally passes 
through the WMAP-favored values, when we go to higher gaugino fractions Zg. For the 
Ghana*^"-* scan (in green) only, the Higgs mass drops below the LEP bound before the 
WMAP-favored region is reached, as is shown by the jagged vertical line. 

8 Summary and outlook 

In this paper, we considered effective field theory corrections to the MSSM, and incorpo- 
rated them in existing packages for accurately calculating the relic density of dark matter. 
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We found that the corrections can make a difference for dark matter in certain regions of 
parameter space, and we performed a first scan of parameter space of these models. 



An important question that we addressed (but did not fully answer) in sections 7.2, 7.3 



and 7.6 is whether the new models generated by turning on the BMSSM corrections ei^2 
can effectively be recreated by modifying the other pre-existing parameters of the ordinary 
MSSM, like mo or fi. To be concrete, take the BMSSM models in the benchmark scans in 
the slice through parameter space in figure |5| Let us call the models in the red BMSSM 
curve that are not on top of some parts of the black MSSM curve "new BMSSM models" . 
Do ordinary MSSM models that are in some sense physically similar to the new BMSSM 
models exist anywhere in some other slice of parameter space than that fixed in figure 
|5]? If they do, then whether they should be considered equally interesting as the BMSSM 
models that satisfied the constraints "naturally" in this slice depends on the details. For 
example, as argued in section 7^, there could be only a restricted subset of parameters 
that could be varied to reach those MSSM models, and those parameters could change e.g. 
accelerator physics of these models in characteristic ways. Or, it could require finetuning 
of parameters so awkward (see e.g. section 7.3) that the price of reaching these models in 
the MSSM is not worth paying. This will of course depend on one's idea of naturalness. 
We leave this question here, and hope that our detailed remarks in previous sections give 
some useful guidance for further studies of the importance of naturalness for MSSM vs. 
BMSSM dark matter. 

There are many directions one could take this further. First, we imposed many restric- 
tions on the models, some of which are well motivated, some of which could be relaxed. 
The most obvious direction is perhaps to perform large dedicated scans to study cer- 
tain characteristics of these models more systematically. It is generally challenging to 
give precise boundaries in parameter space of the ordinary MSSM where the model is 
ruled out by experiment (and thus can be ruled in by the BMSSM), since one can often 
tune MSSM parameters to very special points to evade a given bound. Thus, to really 
understand what BMSSM models are truly "new" and do not exist anywhere in MSSM 
parameter space (especially in models more general than the MSSM- 7), a Markov Chain 
Monte Carlo (MCMC) analysis might be more enlightening than the random grid scans 
we have performed here. The light Higgsino region would be an example that would be 
worth studying further in this respect. 

Next, one can of course generalize our rather "minimal" BMSSM implementation fur- 
ther. For example, it might be interesting to study CP violating processes, allowing for 
imaginary parts of the ei,2, or physics directly related to the operators. One com- 

ment about this last point: as discussed in section 4.3, it is interesting that the 
operators in the Kahler potential are rather stringently restricted by the precision elec- 
troweak observables S and T |37j. In the effective theory, the coefficients of those 
operators [7] are independent of the coefficients ei^2 we have considered, but in an un- 
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derlying theory there may be relations between and (see [Tj for an example). As 
a highly optimistic scenario, consider a future experiment where dark matter has been 
detected that is well described by the neutralino in the light Higgsino region of our plot. 
Let us interpret that as evidence for new physics beyond the MSSM. By the logic in this 
paragraph and section 4^, constraints from electroweak precision observables could be 
used to further constrain the underlying theory in this situation. 

Also, one could study the interpretations of specific experimental/observational data 
in the light of these corrections; some recent studies of accelerator physics that we think 
could be relevant to combine with further work on our dark matter calculations include 
[381 EHl Sni SH 112]. Also the work on BMSSM baryogenesis |37j would be interesting to 
study further in the light of our calculations. 

Finally, especially in BMSSM models with more new parameters, it would be useful 
to try to obtain additional restrictions on those classes of models to keep the analysis 
practical. A related question is if it is worth pushing the size of the BMSSM corrections 
further than we have. Suggestions in this direction include allowing for the creation of 
new vacua, but imposing the BDH criterion |27j, and thereby trying to make sure the 
model is stable at least on timescales of the age of the universe. 



A Notation and conventions 

We denote the two Higgs doublets with hypercharge Y = +1/2 and F = — 1/2 by Hu and 
respectively. In the literature [HI |13] one also often finds the notation H2 and Hi, 
i.e., Hu = H2 and Hfi = Hi. In SU{2) components the doublets are written 

Hu=(ii), H,= (^"). (79) 



All suppressed SU{2) doublet indices are contracted in the main text. Doublets trans- 
forming in conjugate representations are contracted with the Kronecker delta 6ab, for 
example we have 

HiHu = SatiH^ynt = \H:\' + \H'S (80) 

while doublets transforming in the same representation are contracted with the epsilon 
symbol eat, with ei2 = —1, for example we have 

HuHd = €abH'^H^ = H^H^- H:^H^. (81) 

It should be clear which contraction is intended where. When multiple contractions appear 
we use brackets to denote which fields to contract, for example we have 

{HuH,){HiHu) = e^b^cdKH'AHlYHi . (82) 
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Our sign of fi agrees with the sign in DarkSUSY and is opposite to that of DST [7]. The 
parameters ei and €2 that we introduced in (31) and (33) have the same sign as in DST; a 
positive ei correction increases the hghtest Higgs mass rriho, as can be seen in figure (ll]). 



A.l Superpotential interactions 



The dynamical field content of a left-chiral super field $ is a complex scalar and a 
left-handed two-spinor ip. Given a superpotential W = W{^i) for a family of left-chiral 
superfields we calculate the fermionic interactions using 



(fcrmion)-^ 



2 I ^ d^id^i 



and the F-term scalar potential using 



tpitjjj + h.c. 



^3) 



dW 



These formulas hold for any W, not just renormalizable ones. 



B Higher-dimensional operators 

There are two obvious questions about effective field theories with higher-dimensional 
operators: how we compute with them, and where they come from. In this short appendix, 
we try to elucidate this as briefly and as simply as possible. Some of this is textbook 
material [SI Ch. 12.3], some is available in review articles like [121 US]- 

Consider the following "microscopic" toy-model Lagrangian for two neutral scalars 
and H^, and three neutral fermions S, and H^: 

■^microscopic = -^bce + '^frec,S ~ H^^H^ — Xs S H^H^ + h.C. (85) 

where -^free are the usual kinetic and mass terms for H^, and and ^iree,s are 
the kinetic and mass terms for the fermion and we introduced the Yukawa coupling A5, 
that we set to be equal for the two terms. We want to calculate an effective Lagrangian 
=Sfcff that incorporates the lowest-order interactions between and due 

to the exchange of S", but no longer contains S explicitly. In other words, we want 
to approximately integrate out S. (We write "approximately" to distinguish this from 
performing the functional integral over S", which would not be an approximation.) 

To this end, recall how the four-fermion coupling in Fermi theory is a low-energy 
approximation of four-fermion interactions in the electroweak theory, valid for external 
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momenta much smaller than the mass of the exchanged gauge boson. Analogously, the 
theory ( 85 ) gives rise to the following operator in the effective Lagrangian at tree level: 



S 




H° 



On the other hand, as in the body of the paper (see eq. (28)) we can write a general 
effective dimension-5 operator (the little circle in the second diagram above) of this form 
with a dimensionless coupling A: 



^6) 



The 1/M suppression for some mass scale M is just by dimensional analysis. Here, the 
energy scale M is supposed to be the "scale of new physics", the energy at which our 
effective theory starts to deviate from the microscopic theory. For the microscopic theory 
above, clearly this will be M ~ Mg. (Whether M is exactly this scale Mg or somewhere 
close to this scale is a matter of definition — this is equivalent to stating when <^ 
fails to be satisfied to a prescribed accuracy.). We now match the effective field theory 
coupling A to the microscopic coupling A5 at a prescribed renormalization group (RG) 
scale A, say A = Ma: 



A(A) = A|(A) at A 



Ms- 



17) 



The couplings may run differently below and above this scale, since the field content is 
effectively different, so the P functions will be different. We see that it is useful to think 
of A (A) and A5(A) as two different quantities, that are related by A = A| only at the 



matching point in (87). 



In our examples in the paper, we express the coupling in terms of ei := — A/i*/M (see 



eq. (32)), so equivalently 



dim— 5 



HuHdHuHd 



To summarize, the effective theory valid below the energy E = M is 



eff 



free 



dim— 5 



^9) 



from (85) from (|86| 
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where 5* does not appear at all. We distinguish this from the microscopic theory (85) 
where S appears explicitly. If we use to compute the cross section a^HuHd 
we have the low energy expansion 



HuHd 



5) 



E 
M 



02 + . 



(90) 



where ctq is the cross section for ei = 0, a\ is the leading correction, and each term in the 
low-energy expansion itself has a loop expansion (do = o"o,o + Acro,i + •••)■ W^e find that 



for E <C M, the cross section (90) with just the linear term is not much different from 



the complete result in the underlying theory -^microscopic- 

At loop level, as long as we use a mass-independent scheme such as MS, the expansion 
in [E/M) will remain consistent. This does introduce apparent divergences, but they 



cancel when we match the theory to the microscopic Lagrangian (85) at the RG scale 
A = M, by a relation similar to (87) (see e.g. [46]). In the same vein, apparent UV 



divergences due to hard supersymmetry breaking operators like ( 33 ) are no more harmful 



than soft supersymmetry breaking terms when matching to the microscopic theory, where 
supersymmetry breaking is usually spontaneous (one clear example of this in the MSSM 

is m)- 

So to answer the questions posed at the beginning of this appendix, we calculate with 
the higher-dimension operators as with any other Lagrangian, as long as we remember 
that the approximation is only valid at energies much smaller than the scale of new physics 
M, which can be set to in this example. It is also important that in the effective field 
theory, we can work without direct reference to any parameters in -^microscopic, and view 
ei in -^eflf as a model parameter. 

Now for the origin of the operators, and the generality of the effective theory, versus 
that of the microscopic theory. We have seen that for a given microscopic theory, we can 



express the low-energy parameters in terms of microscopic quantities, as in (87). It is 



important, however, that (85) is not the only theory that can give rise to new operators 



that can be approximated by the form (89). For example, if we instead integrate out a 



field with SU{2) charge, we can obtain a similar expression for -Sfes- The effective theory 



captures aspects of several possible underlying theories, not just eq. (85). In this sense it 
is more general than any particular microscopic theory. 

The effective theory is also less general than any particular underlying theory, in the 



following sense. The particular Lagrangian in (85) is actually part of the Next-to-Minimal 



Supersymmetric Standard Model, or NMSSM. In this theory, 5* (sometimes called the 
"singlino") can be the LSP in some part of parameter space. That is clearly not possible 



in the effective theory (89), where the singlino does not appear at all. In other words. 



we only consider indirect effects of any new physics on the interactions between existing 
fields in our theory, we do not consider any new fields beyond the field content of the 
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usual MSSM. In this sense, the effective theory does not capture all aspects of a given 
microscopic theory. 

One final comment on the range of validity of the effective theory. If we require the 
microscopic theory to be perturbative, e.g. A < 1, then if we view ei as fixed at say 0.01, 
we find a constraint 

\a\ ei 0.01 , , 

F = x = ^>°-" 

so for example for M ~ 5 TeV, we have /i > 50 GeV for perturbativity. Similar comments 
hold for €2 and rrisusY- For practical purposes, we will allow M to vary slightly, such 



that (91) is always satisfied for our (reasonably small) range of the parameters and 
"^susY- Thus, we do not bring this up in the discussion in section |4.3[ One can also 
contemplate viewing the effective operators as not arising from a particular perturbative 
microscopic theory at all, as for the chiral Lagrangian in QCD, but we would like to keep 
perturbativity as an option. 



C Feynman rules 



In this appendix we present the Feynman rules generated by the Lagrangian in eq. (37). 
These include interactions between charginos, neutralinos and Higgs bosons as well as 
Higgs boson self-interactions. Concerning the Higgs boson self-interactions, we present 
only three-point couplings and no four-point couplings, as the latter contribute negligibly 
to relic density calculations. 

When a coupling contains two fermions, our convention (following [34J) is that the first 
one always appears with a bar in the Lagrangian while the second one always appears 
without the bar, i.e., the first fermion is outgoing and the second fermion is ingoing. The 
scalar bosons all appear in the couplings exactly as they appear in the Lagrangian. For 
example, to the following Lagrangian term 

^g^H^H.-o^oH^H-^.PLx'j (92) 



we assign the following coupling: 



Here we assumed that the coupling is symmetrized in 

All Feynman rules contain an i that we do not write explicitly, i.e. all couplings below 
should be multiplied by a factor of i. The DarkSUSY convention is to have the vertices 
divided by a factor of i, so that the vertices below are coded into DarkSUSY exactly as 
written. 
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C.l H-H-H vertices 

The vertex is 



H 

/ 



I 

I 

/ 

/ 

H 4 

V 

\ 

\ 

V 

\ 

\ 



9hhh 



H 



9h°h°h° — —^iv ^12-\/2sin(/3 + a) — 6-\/2cos2Q;sm(/3 — a)^ 

+ ^— 3\/2cos(/3 — a) + 3\/2cos2q;cos(/3 + a)^ 
9h°h°h° — —^iv ^6-\/2cos2q;cos(/3 — a)J 

+ (^V2sm{P — a) — 3A/2cos2ci;sin(/3 + a)^ 
9h°h°h° = -^iv (^6V2cos2Q;sin(/? - a)^ 

+ e2V ^v^cos(/3 — a) — 3-\/2cos2q;cos(/9 + a)^ 
9h°h°h° = —^iv ^12\/2cos(/3 + a) — 6\/2 cos2q; cos(/3 — a)^ 

+ e2V (^-3V2sin(/3 - a) + 3V2cos2Q;sin(/3 + aj^ 
9hOhOho = -^iv (^-2^/2 COS 2/3 sm{p - a)^ 

+ e2V (3V2cos{p - a) - V^cos2/3cos(/3 + a)) 
9h^h^go = —eiv ^2-\/2cos2/3cos(/? — a)^ 

+ e2V (y2sin{(5 -a) - V2cos2/3sin(/3 + a)) 
9h^g°g° — —^iv ^4-\/2sin(/3 + a) + 2v^cos2/3sm(/3 — a)^ 

+ e2V (^-v/2cos(/3 - a) + \/2 cos2/3cos(/5 + a)^ 
9h°h°h° = -^iv (^2^2 cos 2/3 cos (/3 - a)) 

+ e2V (^3v^sin(/'j - a) + ^2 cos 2/3 sin(/3 + a)j 
9hoh°go = -^it^ (^2V2cos2/3sin(/3 - a)^ 

+ e2V ^-v^cos(/3 - a) - cos 2/3 cos (/3 + a)) 
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Qh^g^g^ = —eiv ^4\/2cos(/9 + a) — 2^/2 cos 2/3 cos(/9 — a) 

+ 62^ (-\/2s\u{l3 - a) - V2cos2/3cos(/? + a)^ 
Qh^h+h- — —^iv ^— 2 cos 2/3 sin (/3 — a)^ 

+ ^\/2 cos(/3 — a) — \/2 cos 2/3 cos(/3 + a 
9h°h+g- — —^iv ^2-\/2cos2/3cos(/3 — a)^ 

+ e2V (^-v^cos2/3sin(/3 + a)^ 
9h°g+h- — —^iv ^2-\/2cos2/3cos(/3 — a)^ 

+ e2V (^-\/2cos2/3sin(/3 + a)) 
9h°g+g- — —^iv ^4-\/2sin(/3 + a) + 2-\/2cos2/3sin(/3 — a 

+ e2V (^-\^cos(/3-q;) + \^cos2/3cos(/3 + q;)) 
9h°h+h- — —^iv ^2\/2 cos2/3cos(/3 — a)^ 

+ 62^; (^V2sin(/3-Q;) + \/2cos2/3sm(/3 + Q;)) 
9h^h+g- = -^iv (^2v^cos2/3sin(/3 - a)j 

+ 62^ (^-V2cos2/3cos(/3 + a)^ 
9H0G+H- = -eiV (^2A/2cos2/3sin(/3 - a)j 

cos2/3cos(/3 + a)^ 
9h°g+g- — —^iv ^4-\/2cos(/3 + a) — 2-\/2cos2/3cos(/3 — a) 

+ e2V (^-V2sm{(3 - a) - \/2cos2/3sin(/3 + a)) 

9H0H+G- = ^2^^ (-iV2) 

9h0g+h- = ^2V [iV2^ 



C.2 H-x-X vertices 



C.2.1 H°-XrX°i vertices 



The vertex is 
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9ho^o^o = 2 ( ~ vV2 cos P cos aN*^N*4 - vV2 sm Psm aN*^N*^ 

- 2V2v sm{a + P)1{N*,N;, + TVJ^TV* )) (114) 

- 2V2v cos(a + /?)^(iV;^A73 + N;,N:,)) (116) 



4°x?x? = 2 (^-^) ( - ^^71 «i^2/?7Vi^A7, - z^— sin 2/3iVi^iV;3 

- t2V2v^{N:,N;, + N;,N:,)) (118) 

3Him°j ^ fex»X») (119) 

fe=fex")* (120) 



C.2.2 H--^j-xt vertices 

The vertex is 
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^H-X'xt 



//* 

e" 
/i 
ei 



(-2^; cos^ pN*^V*^ - V sin 2/?Ar;3i;* ) 

f sin 2l3Nj^Uc2) 
V sin 2(5N*^V*^ + 2v sin' /3A73K*2) 



- sin 2f3Nj3Uc2 + 2v cos' (3NjiUc2) 



Changing the directions of the arrows gives 



H- -■ 




C.2.3 H^-xt-Xd vertices 
The vertex is 




^HOx+xt-^^ + dHo^tx^PR 





































(-1 



V2vsm{a + P)U:,V:, 



'c2 



V2vcos{a + (3)U:,V:, 
V2v cos{a + P)Ud2Vc2 
45 



^kmi = (-^) ^^^-ut^y^. (129) 

4o^+^+ = {-'^ HV2v)Ud2V,2 (130) 



C.3 H-H-x-x vertices 
C.3.1 H^-H^-xi-x^j vertices 
The vertex is 




-sm2a^{N*,N;, + N;,N:,)) (131) 

+ sm2a^{N:,N;, + NJ^N*,)') (133) 

.W.,^.(-^)(^..«.^o.« 

+ sm2/3^(7V*7V;3 + 7V;^7V*)) (135) 

3HiHix°x', = te^^§x°x°) (^36) 

- sm2p^{N:,N;, + N;,N:,)) (137) 
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-2{-—\(lsin2aN* N* - I sin 2aN*N: 



- 2cos2a^(iV*,iV;3 + iV^.iV^,)) (139) 

9h°h0x°x° ^ ^ ^) ^ — i COS a sin (3N*^N*^ — i sin a cos (3N*^N*^ 

- 2icos{a - (3)\{N;,N*, + N*,N;,)) (141) 

+ 2isin(a - P)\{N:,N*^ + 7V;4A^i^)) 

9hohOxOx° = 2 (~^) (^sinQ;sin/3A^*4Ar*j -icosQ;cos/3A^*3A^*3 

+ 2zsin(a - P)1{N:,N;, + N;,N:,)) (144) 

9H°Hix°X° = {9hOH°x^X°) (14^) 

9H°G°^i^j " ^ (~ ^) ( ~ ^ P^u^ji -icosa sin (5N*^N*^ 

+ 2zcos(a - /3)^(7V;iA73 + ^^4^^;.)) (1^6) 

4°G0x?x? = 2 (- ^) ( - ^ sin 2/?iV;^iV;, + \ sin 2/3iVi;, A73 

- 2 cos 2/i^(i\r;^A73 + iV;4iV^,)) (148) 



C.3.2 //+-i/--x?-x° vertices 

The vertex is 
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H- 



H+ 




9G+G-xm - 2 



9H+H-xm 



9G-H+xm 



9G-H+xm 



SG+H-xm 



9G+H-xm 



/I* 



sm2P-{N*^N;, + N;,N:,) 



9G+G-xm 



- [9H+H-xm 



1^ 



-2cos'P-{n*,n;, + n;,n:,) 

+2 sin' P^iNi^Nja + Nj^Ni^) 
+2sm'p^{N:,N;, + N;,N:,) 
-2 cos' pliNi^Njs + Nj^Nis) 



(150) 
(151) 
(152) 
(153) 

(154) 

(155) 
(156) 
(157) 



C.3.3 H^-H^-xt-Xd vertices 
The vertex is 




9H°H0xtx^ - 2 



sm2aV:,U:, 



(158) 
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dioH^xtxi 



L 



9gogox+x^ 



^HlH^txi 



^HlH^ixtxt 



9H°Hixtxt 
dioH^xtxi 
9H°G0x+xi 
^HfGOxtx^ 



^HiiHUtxi 



^H^Hixtx^ 



L 

dn^GOx+xt 



dnOGOx+xi 



^H^GOx+xi 
^HlG^xtxi 



£1 
//* 

f* 

/i 

/X* 

A 



sui2aVc2Ud2 



1 



- sm2aVc2Ud2 

^ sin 2/3^2 C/d2 
sm2/3V:,U:, 

sm2pV,2Ud2 



fi* 

(=* 
._L 

fi 

ei 
fi* 

A 

//* 

A 
1^ 

/i* 

.1 
/i 

ei 

/X* 

£1 

ei 

._L 
1^ 



cos2aV:^U:^) 



cos 2Q;V;2C/d2) 



^ cos(/3 - a)y;2C/:2) 

-icos(/? - a)Vc2Ud2) 
tsm{P-a)V:,U:,) 
-ism{P - a)Vc2Ud2) 
tsm{P-a)V:,U:,) 
-isin(/3 - a)Vc2Ud2) 
-icos{P-a)V:,U:,) 
+i cos(/3 - a)Vc2Ud2) 



cos2pV^,U:,) 



cos2/5K2^d2) 



(159) 
(160) 
(161) 
(162) 
(163) 
(164) 
(165) 
(166) 
(167) 
(168) 
(169) 
(170) 
(171) 
(172) 
(173) 
(174) 
(175) 
(176) 
(177) 
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C.3.4 H+-H--xt-X^ vertices 
The vertex is 

H- 



H+ 




9h+h- 



L 

9g+g- 



9g+g- 



Xc Xd 



Xc Xd 



Xc Xd 



G-H+xtxJ 

B. 

G-H+xtxt 



9g+ 



H-Xtxt 



-,R 



9G+H-xtxi 

C.3.5 H--H--{xtT-Xd vertices 
The vertex is 



H- 




/i* 



/i 
/I* 



e 



{-2sm2pv:,U:,) 

(-2 sin 2/3^2^/^2) 
{2sm2pv:,U:,) 

{2sm2i3VM2) 



M (-4sin2/3K2C/d2) 



(-4 sin^/3K* 



(178) 
(179) 
(180) 
(181) 
(182) 
(183) 
(184) 
(185) 
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9h- 


H 


p 




9h- 


H 


T 




Sa- 


G 






ul 


G 








H 








H 



4 - 



4 - 



Changing the directions of the arrows gives 

H- ixtr 



II* 

/i* 

.!i 



— sm 



cos 



-sin2/3C/c2C/d2) 
sin2/3K2Kr2) 

-Sin2/3t/e2t/d2) 




C.3.6 H^-H--x°-xt vertices 
The vertex is 



H- 




^HiG-mi 



9H0G-mi 



(~^) (-^cosacos/?^;^^*^ - v^cos/^sinaAT^gl^^^) 
— — j ^\/2 sin a sin (3Nj^Uc2 + cos a sin (3Nj4JJc2^ 
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-"2-" XjXc Y n 



V2 COS a sin pN*^V*2 + \/2 sin a sin /^A^gl^^a) 
g^Q^_~o~+ = (^~~^ (y2smacosPNj3Uc2 + V2cosacosPNjiUc2^ 
silG-x'.xt = {~^) (\/2sinQ;cos/3Ag^K2- v^cosacos/^A^aKa) 
dnOG-x'^xt ^ (^^°^'^^™'^-^j3C42 - V^sinQ;sin/3A^j4C/c2) 

(^- V^sin a sin PN*^V*2 + ^2 cos a sin 
\/2 cos a cos (5Nj^Uc2 — sin a cos (5Nj/JJc2j 
9H0G-~^~xt = (-^) (-^^sin2/5iV;,V;* -zV2cos^/?iV;3K* 
4oG-3^o^+ = (^-^) (^-^-^sin2/3A,-3t/e2-iV2sin2/3iV,-4^7e: 
^Vx?x^ = (-^) (^^sin^/37V;,K2 + ^^sm2/3iV;3K2) 
^Sfi^-x^x^ ^ (-^l) (-^^cos'W-3C/c2-?^sin2/3A,-4C/c2 
4o^_^o^, = (-^) (^+iV2cos^/?7Vj4K2-i^sin2/?7V;3K*, 
9aoG-x!^.^+ = (^-^) (^-i\/2sin2/3Ar^-3t/e2 + i-^sin2/3Ar^-4t/c2^ 
5SoH-^o,+ = (-^) (-^-^sin2/3iV;,K2 + ^V2sin2/3iV;3V;*2^ 
^2o^_^o^+ = (^-^^ (^-i-^sin2/3iV,3[/,2 + iV2cos2/3A^,-4C/c2 



''c2 



Changing the directions of the arrows gives 
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